Information transfer in stochastic optimal control theory with information theoretic criterial and application

ABSTRACT

The current disclosure provides a method for transmitting encoded information signals through a control system and to a decoder. The encoded information signals are transmitted along with control signals as an encoded message. The information signals are encoded based at least in part on a control-coding capacity of the control system.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation of U.S. patent application Ser. No. 15/461,462, which was entitled “The Information Transfer in Stochastic Optimal Control Theory with Information Theoretic Criterial and Application” and filed on Mar. 16, 2017 which claims priority to and the benefit of the filing date of U.S. Provisional Patent Application Ser. No. 62/309,315, which was entitled “The Information Transfer in Stochastic Optimal Control Theory with Information Theoretic Criterial and Application” and filed on Mar. 16, 2016. The entire disclosure of this application is hereby expressly incorporated by reference herein for all uses and purposes.

TECHNICAL FIELD

This patent application generally relates to coding of information signals in controllers of general control systems, and using the control system to communicate the encoded information to processes connected to control systems, such as field devices, outputs, other controllers or processors. The invention further relates to: (i) controller design in control systems with a dual role, (ii) achieving stability and optimal performance with respect to control objectives, (iii) encoding information, (iv) transferring the information through the control system, and (v) decoding or estimating the information at outputs of control systems, with arbitrary small error probability. The invention facilitates communication from one processor of a control system to another processor.

BACKGROUND

Current telecommunication and information technology systems are designed based on Shannons operational definitions of coding-capacity, the maximum rate of communicating information over noisy channels, and coding-compression, the minimum rate of compressing information, which give the fundamental performance limitations for reliable communication. They utilize encoders and decoders, to combat communication noise and to remove redundancy in data.

Current dynamical control systems are designed by utilizing feedback controllers, actuators and sensors, to ensure stability, robustness, and optimal control objectives and performance.

Industries related to modern control and communication systems have experienced tremendous growth due to their increasing applications in information technology which affects everyday lives of people. The next generation of engineering systems integrates control, communication, protocols, etc. to develop complex engineering systems, which can be implemented in energy systems, transportation systems, medical systems, surveillance networks, financial instruments etc. Many of these applications consist of multiple control sub-systems and communication sub-systems integrated together to achieve control and communication objectives.

In the field of communication most systems are designed by utilizing encoders and decoders, to combat communication noise and to remove redundancy in data. Current telecommunication systems, whether point-to-point, network, mobile, etc., are designed, based on Shannon's operational definitions, which give the maximum rate of communicating information over noisy channels, and the minimum rate of compressing data generated by sources, called coding-capacity of communication channels, and coding-compression of information processes, respectively.

In the field of modern control theory and applications, controllers are designed to control the control system outputs, to optimize performance and to achieve robustness with respect to uncertainties and noise. Over the year the criterion of optimality of controllers is the average of a real-valued sample path pay-off functional, of the control, state and output processes, which is optimized to achieve optimal control performance. It is well-known that controllers are designed to control the control systems but not to encode information and communicate this information, through the control system to other processes.

The general separation of control system design and communication system design has divided the community of developers into independent groups developing controllers for control systems, and communication systems for noisy channels. Shannon's operational definitions of achievable coding-capacity, and coding-compression, which utilizes encoders and decoders, to combat communication noise and to remove redundancy in data, did not apply beyond communication system, such as, in the ability of control system to transmit information.

To this date, no controllers in a control system is designed to communicate information, from one controller, through the control system, to another controller, or to any of its outputs, or from elements of a control system to elements of another control system, using encoders and decoders as done in problems of information transmission over noisy channels.

To this date, Shannon's operational definitions of coding-capacity and coding-compression are not used beyond communication system, such as, in dynamical control system. It is well-known that encoders are designed to encode and transmit information over noisy communication channels but not to control channel output processes.

SUMMARY

In an embodiment, a method of designing an controller-encoder pair in control systems or general dynamical systems or decision systems, comprising the controller controls the control system, the encoder encodes information signals and transmits the information over the control system to any element attached to it, and reconstruct the encoded signals using decoders, with arbitrary small error probability.

The method described above further includes determining the expression of Control-Coding Capacity (CC-Capacity) of the control system, for stable or unstable systems, wherein the expression represents the maximum rate in bits/second of simultaneously controlling the control system and encoding information, and decoding the information at any processors attached to the control system, by solving the expression of CC-Capacity.

The method of above further includes simultaneous

(i) optimal performance with respect to control objectives, and

(ii) optimal coding of information signals with respect to communication objectives.

The method of above further includes the the specification of minimum power necessary to meet the optimal performance with respect to control objectives, and the specification of excess power which is converted into transmitting information through control systems with positive data rates.

In another embodiment, the controller-encoder comprise of deterministic part and a random part, wherein the deterministic part controls the control system, while the random part encodes information.

The method further includes determining how existing controllers which are not designed to encode information, can be modified to encode information, and transmit information over control systems, and decoded at any processors attached to the control system, such as, in control-to-control communication, etc.

BRIEF DESCRIPTION OF DRAWINGS

FIG. 1 is a block diagram of an example of controller-encoder of a control system, which simultaneously controls the control system and encodes an information signals, and transmits the information signals over the control system, and decodes the information signals at the output of a decoder and/or estimator;

FIG. 2 is a flow chart depicting a method of an example method for transmitting encoded information signals and control signals in a control system.

DETAILED DESCRIPTION

The techniques of the present disclosure facilitate a method of controller design in control systems or decision systems, which is based on utilizing controllers, encoders and decoders, wherein the encoders and decoders are designed using techniques from communication system designs. Further, these techniques may be universal to any dynamical system with control inputs and outputs. The maximum rate for simultaneous control of the dynamical system and transmission of information over the dynamical system is called the Control-Coding Capacity (CC-Capacity).

The design of the controller-encoder for control systems, aims at simultaneous control of the control system and communication of information from one processor to another processor. More specifically, the CC-Capacity of dynamical control systems identifies the interaction of control and information, and controller-encoders-decoders are designed to achieve simultaneously optimal performance with respect to control objectives, and optimal transmission of information, via the dynamical system, from one controller to another controller or output. etc.

Applications of the CC-Capacity to control systems, and in general, to any dynamical system with inputs and outputs, include biological, financial, quantum, etc., because any such system has CC-Capacity.

To this end, a methodology determines how to design control systems, which in addition to effectively control system outputs (e.g., stabilize them), they also encode information, and communicate the information from one controller to any other controller, via the control system which acts as a communication channel, and the information is decoded or estimated at other processors attached to the control system, with arbitrary precision.

Dynamical Systems or Control Systems with Communication Capabilities

FIG. 1 is a block diagram of an example dynamical system or control system 100 including at least some components that may be configured according to and/or utilize the methods discussed herein. In this example system, an information signal 101 is received (e.g., including any suitable information). A randomized strategy also called controller-encoder 102 is designed to control the control system 106 and to encode the information signals 101, the control process or encoded message 105, including a control signal and an encoded information signal, is transmitted as a signal over the dynamical system also called control system 106. A decoder or estimator 108 decodes the process output 107, which is the output of the encoded message 105 after going through the control system 106, to produce an estimate of the information signal 101. The process output 107 is also returned to the controller-encoder 102 as feedback. In some embodiments the process output may be delayed before being returned to the controller/encoder.

The information signal 101 includes analog or digital data. For example, the information signal may include suitable computer generating messages, including digital data, such as digital data representing pictures, text, audio, videos, or signals generated by other control systems to etc. Generally, the information signal 101 may include any number of devices or components of devices generating messages to be encoded by the randomized strategy also called controller-encoder 102 and transmitted over the control system 106. The information signal 101 may include information sources and/or tracking signals as well as any other pertinent data that may be encoded and transmitted through the control system.

The randomized strategy called controller-encoder 102 encodes messages received from the information signals 101 {X₀, X₁, . . . , X_(n)} and uses feedback from the control system outputs {Y₀, . . . , Y_(n)} to generated the control process {A₀, . . . , A_(n)}. The controller/encoder 102 may comprise a control module 103 and a coding module 104. The control module 103 may be used for typical control system activities, such as generating control signals for elements of the control system 106. The control module may also receive feedback 107 and other information related to the control system which may be implemented in controlling the control system. The coding module 104, which is not typical of previous known controller devices, may be implemented in encoding information signals 101 and transmitting the control process 105 as an encoded signal along with control signals to the control system 106. The controller/encoder 102 may implement the control module 103 and coding module 104 in any suitable combination to produce the control process 105. In some embodiments, the control process 105 can be represented as a randomized control signal or a randomized control action.

The coding module 104 may implement any of the techniques described below to encode the information signals 101. For example, the encoding module may be designed based on the capacity of the control system, known as the control-coding capacity, of one or more control systems 106 and then encode the signals based on the determined capacity. Further, the coding module 104 may also implement any known or future developed schemes for encoding the information signal 101.

The dynamical system also called control system 106 may include any number of systems, described by nonlinear recursive models driven by a control process and the noise process, such as, the dynamics of a ship, submarine, missile, plane, car, and any model of a dynamical system which is controlled by a controller, such as a process control system. The output of the control system 107 is available to the decoder or estimator 108 and comprises encoded information signals and control signal feedback. The decoder 108 may be any device suitable for receiving the encoded information signals from the process output 107 such as a sensor, field device, another controller, etc. In the example of FIG. 1 the control system or dynamical system 106 is described by a sequence of conditional distributions {P_(Y) _(i) _(|) _(i-1) _(,A) _(i) :i=0, . . . , n} called conditional distributions or control system distributions.

FIG. 2 is a flow chart of an example method for transmitting encoded signals through a control system. The method starts at block 210 when an information signal is received by a controller of a control system. The information signal can be of any form discussed above with regard to FIG. 1. At block 220, a controller may determine one or more control signals that to be sent to the control system. At block 230, the received information signals may be encoded. The information signals may be encoded by the controller at a coding module, for example. The information signals may also be encoded based on any of the techniques discussed herein. Next, at block 240 the controller may transmit the encoded information signals along with the control signals through the control system.

Previous systems would only transmit control signals through control systems. These previous systems lacked the capability of sending information signals in such a manner. For example, if information needed to be transmitted through a control system or from one controller to another controller in the control system, this was not enabled by the previous systems. In-stead, previous systems would require additional communication channels for communicating information signals and they do not use the control systems as communication channels used to transmitting information signals. In other words, by transmitting encoded messages which include both control signals and encoded information signals through the control system, the current system enables simultaneous control and communication of data transmission. Further, the current system requires further utilizes the control system as a communication channel and also provides more optimized communications through currently available control systems. The implementation of such techniques was not previously considered (using communication channel optimization techniques in a control system).

The encoded message transmitted in block 240 may be utilized at two different junctures because the encoded message essentially contains two parts, an encoded information signal and one or more control signals. At block 250 the encoded information signal portion may be decoded. The encoded portion may be decoded by a decoder as described above with respect to FIG. 1. In an embodiment, the encoded signal may be decoded using a particular preset encoder scheme. In another embodiment, the encoded signal may indicate the encoder scheme to be used for decoding the message.

At block 260, the control signal portion may be received by the control system. The control signal portion may be used to control one or more elements of the control system. In an embodiment, once the control signal is received, feedback may be generated by the device or devices that have received and utilized the control signal. For example, feedback may be generated and transmitted as part of process output, such as process output 107 discussed above with respect to FIG. 1. The control signals may be received by a controller, a field device, a sensor, or any other suitable device in a control system.

Although blocks 250 and 260 are discussed in sequential order, these steps may take simultaneously or in any suitable order. In an embodiment, the present techniques allow the control system to be used as a communication channel for transmitting an encoded information signal from a controller to another device which may decode the encoded information signal.

All elements of dynamical systems or control systems 106 with communication capabilities shown in FIG. 1 are analogous to elements source, encoder, channel, decoder of transmitting information through noisy communication channels, as explained below.

-   -   (i) A^(n) Δ{A₀, . . . , A_(n)} is the control process of the         control system, and A^(n)=a^(n) ∈x_(i=0) ^(n)         _(i) ^(n) the control actions. This is analogous to channel         inputs in problems of information transmission over noisy         communication channels.     -   (ii) Y^(n) Δ{Y⁻¹, Y₀, . . . , Y_(n)} is the output process of         the control system taking values in x_(i=−∞) ^(n)         _(i), and Y⁻¹ is the initial state. This analogous to channel         outputs in problems of information transmission over noisy         communication channels.     -   (iii) {P_(Y) _(i) _(|Y) _(i) _(-1,A) _(i) :i=0, . . . m} is the         control system conditional distribution. This is analogous to         the channel distribution of the communication channel in         problems of information transmission over noisy communication         channels.     -   (iv) {         _(A) _(i) _(|A) _(i-1) _(,Y) _(i-1) : i=0, . . . , n} is the         randomized control strategy or conditional distribution of the         control process. This is analogous to channel input conditional         distribution in problems of information transmission over noisy         communication channels.     -   (v) X^(n) Δ{X_(i): i=0, . . . , n} is the information source or         tracking signal taking values in x_(i=0) ^(n)         _(i). This is analogous to information source and X^(n)=x^(n)         the specific message in problems of information transmission         over noisy communication channels.

In view of (i)-(v) there is a one-to-one analogy between elements of control systems and elements of noisy communication channels in problems of information transmission.

A. Consider a general control system with control process A^(n) Δ{A_(i): i=0, . . . , n}, output process Y^(n) Δ{Y: i= . . . , −1, 0, 1, . . . , n}, control system conditional distribution P _(Y) _(i) _(|Y) _(i-1) _(,A) _(i) ≡Q _(i)(dy _(i) |y ^(i-1) ,a ^(i)),i=0, . . . ,n

a set of randomized control strategies or conditional distribution of the control process

_([0,n]) Δ{P _(A) _(i) _(|A) _(i-1) _(,Y) _(i-1) ≡P _(i)(da _(i) |a ^(i-1) ,y ^(i-1)):i=0,1, . . . ,m}

and a set of randomized control strategies satisfying power constraints defined by the set

${\mathcal{P}_{\lbrack{0,n}\rbrack}(\kappa)}\overset{\Delta}{=}{\left\{ {{P_{i}\left( {\left. {da}_{i} \middle| a^{i - 1} \right.,y^{i - 1}} \right)},{i = 0},\ldots\mspace{14mu},{{n\text{:}\frac{1}{n + 1}{E\left( {\ell_{0,n}\left( {A^{n},Y^{n}} \right)} \right)}} \leq \kappa}} \right\} \Subset \mathcal{P}_{\lbrack{0,n}\rbrack}}$

where κ∈[0, ∞] is the total power with respect to the cost function l_(0,n)(a^(n), y^(n-1))∈[0, ∞).

B. The pay-off of model of A. is the so-called directed information from A^(n) to Y₀ ^(n) conditioned on Y⁻¹, denoted by I(A^(n)→Y^(n)) and defined by

${I\left( {A^{n}->Y^{n}} \right)}\overset{\Delta}{=}{\sum\limits_{i = 0}^{n}{E\left\{ {\log\left( {\frac{{dP}_{{Y_{i}|Y^{i - 1}},A^{i}}\left( {\cdot \left. {Y^{i - 1},A^{i}} \right)} \right.}{{dP}_{Y_{i}|Y^{i - 1}}\left( {\cdot \left. Y^{i - 1} \right)} \right.}\left( Y_{i} \right)} \right)} \right\}}}$

where {P_(Y) _(i) _(|Y) _(i-1) :i=0, . . . , n} is distribution generated from {Q_(i)(dy_(i)|y^(i-1), a^(i)), P_(i)(da_(i)|a^(i-1), y^(i-1)): i=0, 1, . . . , n}.

The Finite Time Horizon (FTH) information CC-Capacity of the control system is to determine an optimal randomized control strategy {P*_(i)(⋅|⋅): i=0, . . . , n}∈

_([0,n])(κ) which maximizes the directed information pay-off defined by

${J_{A^{n}\longrightarrow Y^{n}}\left( {P^{*},\kappa} \right)}\overset{\Delta}{=}{\sup\limits_{\mathcal{P}_{\lbrack{0,n}\rbrack}{(\kappa)}}{{I\left( {A^{n}->Y^{n}} \right)}.}}$

The Control-Coding Capacity (CC-Capacity) called the supremum of all transmission rates over the model of A. is the per unit time limiting version

${C(\kappa)}\overset{\Delta}{=}{\lim\limits_{n\longrightarrow\infty}{\frac{1}{n + 1}{J_{A^{n}->Y^{n}}\left( {P^{*},\kappa} \right)}}}$

C. The method of B. further states that CC-Capacity is analogous to the definition of coding-capacity or capacity in problems of information transmission over noisy communication channels, as follows.

(1) The operational definition of controller-encoder-decoder strategies is a generalization of Shannon's operational definition of coding-capacity in problems of information transmission over noisy channels.

(2) In digital communication theory, the source coding problem can be lossless or lossy (based on quantization), while the channel coding problem is defined with respect to the quantized or compressed representation X^((n)) of the information process {X_(i): i=0, . . . , }.

Operational Control-Coding Capacity of Control Systems or Dynamical Systems

Consider a filtered probability space (Ω,

, {

_(i): i=0, 1, . . . , N},

) on which the following processes and RVs are defined. X ^((n)): ∩

^((n)) Δ{1, . . . ,M ^((n)) },A _(i): Ω

_(i) ,Y _(i): Ω

_(i) ,i=0, . . . ,n, Y ⁻¹: Ω

⁻¹ ,{circumflex over (X)} ^((n)):Ω

^((n)).

A controller-encoder-decoder for a dynamical system, such as, a control system, with power constraint, over the time horizon {0, 1, . . . , n} is denoted by

$G\text{-}{SCM}\text{:}\mspace{14mu}\left( {\mathcal{M}^{(n)},\left\{ {\mathbb{A}}_{i} \right\}_{i = 0}^{n},\left\{ {\mathbb{Y}}_{i} \right\}_{i = {- \infty}}^{n},\left\{ {{P\left( {X^{(n)} = x^{(n)}} \right)} = {\frac{1}{M^{(n)}}\text{:}\mspace{14mu} x^{(n)}\epsilon\;\mathcal{M}^{(n)}}} \right\},\left\{ {g_{i}\left( {x^{(n)},a^{i - 1},y^{i - 1}} \right)} \right\}_{i = 0}^{n},\left\{ {Q_{i}\left\{ {\left. {dy}_{i} \middle| y^{i - 1} \right.,a^{i}} \right)} \right\}_{i = 0}^{n},{\ell_{0,n}\left( {a^{n},y^{n}} \right)},{d_{n}\left( y^{n} \right)}} \right)$

and consists of the following elements.

(a) A set of uniformly distributed messages X^((n)) taking values in

^((n)) Δ{1, . . . , M^((n))}, known to both the encoder and decoder (the controller does not need to know these).

(b) A set of controller-encoder strategies mapping messages and feedback control information into control actions defined by ε_([0,n]) ^(S)

{g_(i):

^((n))×

^(i-1)×

^(i-1)

_(i) :i=0, . . . ,n:a ₀ =g ₀(x ^((n)) ,y ⁻¹),a ₁ =g ₁(x ^((n)) ,y ⁻¹ ,a ₀ ,y ₀), . . . ,a _(n) =g _(n)(x ^((n)) a ^(n-1) ,y ^(n-1)),x ^((n))∈

^((n))}

The set of admissible controller-encoder strategies subject to power constraint κ is defined by

${\mathcal{E}_{\lbrack{0,n}\rbrack}^{S}(\kappa)}\overset{\Delta}{=}{\left\{ {{g_{i}\left( {x^{(n)},a^{i - 1},y^{i - 1}} \right)},{i = 0},\ldots\mspace{14mu},{{n\text{:}\mspace{14mu}\frac{1}{n + 1}{E^{g}\left( {\ell_{0,n}\left( {A^{n},Y^{n}} \right)} \right)}} \leq \kappa}} \right\} \Subset \mathcal{E}_{\lbrack{0,n}\rbrack}^{S}}$ where κ∈[0, ∞) is the total cost or power of the controller-encoder. For any message x^((n))∈

^((n)) and feedback information u _(x)(n)=(g ₀(x ^((n)) ,y ⁻¹),g ₁(x ^((n)) ,a ⁻¹ ,y ⁻¹ ,a ₀ ,y ₀), . . . ,g _(n)(x ^((n)) ,a ^(n-1) ,y ^(n-1)))∈ε_([0,n]) S(κ)

is the controller-encoder strategy. The control-coding book of the controller-encoder strategy for the message set

^((n)) is C^((n))=(u₁, u₂, . . . , u_(M) _((n)) ).

(c) A decoder measurable mapping d_(n):

^(n)

^((n)), {circumflex over (X)}^((n)) Δd_(n)(Y^(n)) such that the average probability of decoding error is given by ¹ ¹ The superscript on expectation. i.e. P^(g) indicates the dependence of the distribution on the encoding strategies.

P error ( u ) ⁢ = Δ ⁢ P g ⁢ { d n ⁡ ( Y n ) ≠ X ( n ) } = 1 M ( n ) ⁢ ∑ ⁢ P ⁢ { d n ⁡ ( Y n ) ≠ x ( n ) | X ( n ) = x ( n ) } ≤ n , ϵ n ⁢ ⁢ ϵ [ 0 , 1 )

(d) Conditional independence holds. P _(Y) _(i) _(|Y) _(i-1) _(,A) _(i) _(,X) _(k) =Q _(i)(dy _(i) |y ^(i-1) ,a _(i)),i=0,1, . . . ,n,∀k∈{0,1, . . . }.

(e) The initial data Y⁻¹ ∈

⁻¹ and the distribution μ(dy⁻¹) may be known to the controller-encoder and decoder.

(f) The encoder-controller-decoder for the G-SCM is denoted by (n+1,

^((n)), ϵ_(n), κ).

The control-coding rate is defined by

$R^{{(n)}\;}\overset{\Delta}{=}{\frac{1}{n + 1}\log\;{M^{(n)}.}}$

(g) Another class of controller-encoder strategies, which are nonanticipative with respect to the information process {X₀, X₁, . . . } is defined by

${\mathcal{E}_{\lbrack{0,n}\rbrack}(\kappa)}\overset{\Delta}{=}\left\{ \left\{ {{\left( {{{{e_{i}\left( {\cdot {,{\cdot {, \cdot}}}} \right)}\text{:}\mspace{14mu} i} = 0},1,\ldots\mspace{14mu},n} \right\}{\epsilon\mathcal{E}}_{\lbrack{0,n}\rbrack}\text{:}\mspace{14mu}\frac{1}{n + 1}{E^{e}\left( {c_{0,n}\left( {A^{n},Y^{n}} \right)} \right)}} \leq \kappa} \right\} \right.$

where ε_([0,n])

{e_(i):

^(i)×

^(i-1)×

^(i-1)

_(i) :i=0, . . . ,n: a ₀ =e ₀(x ₀ ,y ⁻¹),a ₁ =g ₁(x ₀ ,x ₁ ,y ⁻¹ ,a ₀ ,y ₀), . . . ,a _(n) =g _(n)(x ^(n) ,a ^(n-1) ,y ^(n-1))}

Note the following.

(3) By the above definition, the message X^((n)) is taken to be uniformly distributed over

^((n)), hence its entropy is H(X^((n)))=log M^((n)). Moreover, unlike most treatments of capacity of communication channels, which often do not impose cost constraints (except for Gaussian channels in which only power constraint on {A_(i): i=0, . . . , n} are imposed), the control system is not required to be stable, and there is a cost constraint on both {(A_(i), Y_(i)): i=0, . . . , n}.

(4) If the entropy H(X^((n))) is below the capacity of the control system, calculated over the time horizon {0, 1, . . . , n}, then for large enough n, we expect to achieve the control objective of meeting the average power constraint, and the coding objective to reconstruct a randomly chosen message X^((n))=x^((n)) at the output of the decoder with small probability of error. Next, we give the precise definition of an achievable control-coding rate and capacity of the control system.

Control-Coding Capacity of Dynamical Systems or Control Systems

(h) A control-coding rate R>0 is said to be an achievable rate (under power constraint κ), if there exists a sequence of controller-encoder-decoder strategies {(n+1,

^((n)), ϵ_(n), κ): n=0, 1, . . . } such that the random processes (A^(n), Y^(n))≡(A^(g,n), Y^(g,n)) depend on the message X^((n)), and satisfy

${\lim\limits_{n\longrightarrow\infty}\epsilon_{u}} = {{0\mspace{14mu}{and}\mspace{14mu}\underset{n\longrightarrow\infty}{\lim\;\inf}\frac{1}{n + 1}\log\; M^{(n)}} \geq {R.}}$

(i) The operational control-coding capacity of G-SCM under power constraint κ is the supremum of all achievable control-coding rates, i.e., it is defined by C(κ)

sup{R: R is achievable}.

(j) Under appropriate conditions, the CC-Capacity is given by C(K) defined by ( ).

We note the following.

(5) In information theory the operational capacity is often called the coding-capacity of noisy communication channels. Since the objective is to control the control system, in addition to encode information, the operational capacity is called the control-coding capacity.

(6) The definition of CC-Capacity states that if a control-coding rate R is achievable under the power constraint and the control system is operated for sufficiently large n, then we can control the control process {Y_(i): i=0, 1, . . . , n} and reconstruct M^((n))=┌e^((n+1)R)┐ messages at the control system output using the decoder or estimator, with arbitrary small probability of error, for sufficiently large n.

(7) As pointed out below, and in view of the general cost constraint, the cost of control is κ^(0,n) (0) while the cost of communication is κ^(0,n)(C)−κ_(0,n)(0).

D. The method of B. further states the optimal randomized strategy {P*_(i)(⋅|⋅): i=0, . . . , n}∈

_([0,n])(κ) is analogous to the optimal channel input conditional distribution in problems of information transmission over noisy communication channels with feedback encoding.

E. An example of the control system is a General-Recursive Control Model (G-RCM) defined as follows.

$G\text{-}{RCM}\text{:}\mspace{14mu}\left\{ \begin{matrix} {{Y_{i} = {h_{i}\left( {Y^{i - 1},A^{i},V_{i}} \right)}},{0 = 1},\ldots\mspace{14mu},n,} \\ {{Y_{0} = {h_{0}\left( {Y^{- 1},A_{0},V_{0}} \right)}},{Y^{- 1} = y^{- 1}},} \\ {{\frac{1}{n + 1}{\sum\limits_{i = 0}^{n}{E\left\{ {\gamma_{i}\left( {{T^{i}A^{n}},{T_{i}Y^{n}}} \right)} \right\}}}} \leq \kappa} \end{matrix} \right.$

where

Assumption A.(i). {V_(i): i=0, . . . , n} is the noise process independent of Y⁻¹, h_(i):

⁻¹×

^(i)×

_(i) and h_(i)(⋅,⋅,⋅), γ_(i)(⋅, ⋅) are specific functions such as linear or nonlinear, and autoregressive models, etc., with Multiple Inputs and Multiple Outputs, and (T^(i)a^(n), T^(i)y^(n))

γ_(i)(T^(i)a^(n), T^(i)y^(n)), where T^(i)a^(n)⊆{a₀, a₁, . . . , a_(i)}, and T^(i)y^(n)⊆{y₀, y₁, . . . , y_(i)} for i=0, . . . , n is any quadratic or nonlinear function;

Assumption A.(ii). The noise process {V_(i): i=0, . . . , n} satisfies P _(V) _(i) _(|V) _(i-1) _(,A) _(i) _(,Y) ₋₁ (dv _(i) |v ^(i-1) ,a ^(i) ,y ⁻¹)=P _(V) _(i) (dv _(i)),i=0, . . . ,n.

Hence, a G-RCM induces a sequence of control system conditional distributions, and these are not necessarily Gaussian.

F. The analogy between feedback capacity of noisy communication channels and stochastic optimal control problems, with directed information pay-off, introduced in A.-D., and depicted in FIG. 1, states that randomized strategies called randomized control strategies can be used to control the control system and to encode information signals or messages, and communicate the information signals over the control system, which acts as a communication channel, to any of the other processes attached to the control system, and reconstruct the information signals by using a decoder or estimator.

Hierarchical Decomposition: Cost of Control and Communication For any finite n, C_(0,n)(κ)

J_(A) _(n) _(→Y) _(n) (P*,κ) is a concave non-decreasing in κ∈[0,∞), and the inverse function of C_(0,n)(κ) denoted by κ_(0,n)(C) is a convex non-decreasing function of C∈[0, ∞). This implies

⁢κ 0 , n ⁡ ( C ) ⁢ = Δ ⁢ ⁢ inf A i ∣ A i - 1 , y i - 1 , i = 0 , … ⁢ , n : 1 n + 1 ⁢ ( a n → Y n ) ≥ C ⁢ E ⁢ { ℓ 0 , n ⁢ ( A n , Y n ) } ≥ ⁢ inf P ⁢ ⁢ E ⁢ { ℓ 0 , n ⁡ ( A n , Y n ) } ≡ κ 0 , n ⁡ ( 0 ) .

Clearly, κ_(0,n) (0)≡κ_(min) is the optimal pay-off or cost of the control system without communication. Then κ_(0,n)(C)−κ_(0,n)(0) is the additional cost incurred if the control system operates at an information rate of at least C.

The cost of communication is given by

${{\kappa(C)} - {\kappa(0)}}\overset{\Delta}{=}{{\lim\limits_{n\rightarrow\infty}{\frac{1}{n + 1}{\kappa_{0,n}(C)}}} - {\lim\limits_{n\rightarrow\infty}{\frac{1}{n + 1}{{\kappa_{0,n}(0)}.}}}}$

In general, κ(C)−κ(0)>0. In view of the above connections, κ_(0,n)(κ) decomposes into two sub-problems, the optimal control sub-problem and the optimal communication sub-problem, which imposes a natural hierarchical decomposition on any optimization problem C_(0,n)(κ).

In addition, we show that classical stochastic optimal control problems, defined by

⁢J 0 , n ⁡ ( P * ) = inf ⁢ E ⁢ { ℓ 0 , n ⁡ ( A n , Y n ) } . are degenerate optimization problems of the FTH information CC-Capacity ( ).

There are several hidden aspects of the FTH information CC-Capacity J_(A) _(n) _(→Y) _(n) (P*, κ), which include the following.

-   -   i) Randomized Versus Deterministic Strategies. Randomized         strategies         _([0,n]) incur a higher value for the pay-off I(A^(n)→Y^(n))         than deterministic strategies,         _([0,n]) ^(D), that is, when {P_(i)(da_(i)|a^(i-1), y^(i-1)):         i=0, . . . , n}∈         _([0,n]) are delta measures concentrated at         a _(j) ^(g) =g _(j)(y ^(g,-1) ,y ₀ ^(g) , . . . ,y _(j-1) ^(g)         ,a ₀ ^(g) ,a ₁ ^(g) , . . . ,a _(j-1) ^(g))         -   for j=0, . . . , n. Indeed, if randomized control strategies             _([0,n]) are replaced by deterministic strategies then             I(A ^(n) →Y ^(n))=0,∀{P _(i)⋅|⋅):i=0, . . . n}∈             _([0,n]) ^(D).         -   This implies that that for any directed information pay-off,             and hence for the FTH information CC-Capacity, randomized             control strategies are responsible for encoding of             information signals, into the the control process {A_(i):             i=0, . . . , n}, which cosists of a deterministic control             strategy part and a random encoder strategy part. This is             fundamentally different from classical stochastic optimal             control problems J_(0,n)(P*), in which the performance over             randomized and deterministic control strategies is the same.     -   ii) Duality Relations. The inverse of the function C_(0,n)(κ)         J_(A) _(n) _(→Y) _(n) (P*, κ) denoted by is given by the         following optimization problem.

Dual Extremum Problem 1.

⁢κ 0 , n ⁡ ( C ) = ⁢ inf ⁢ { E ⁡ ( ℓ 0 , n ⁡ ( A n , Y n ) ) } ≥ ⁢ κ 0 , n ⁡ ( 0 ) ≡ J 0 , n ⁡ ( P * ) .

The inequality states that it costs more to simultaneously control and transmit information than to control only. The additional cost to communicate is κ_(0,n)(C)−κ_(0,n)(0); this is quantified in the application examples.

Moreover, if the randomized control strategies

_([0,n]) are restricted to deterministic strategies

_([0,n]) ^(D), then and necessarily C=0. The resulting optimization problem reduces to the following classical stochastic optimal control problem (without an information theoretic constraint).

Degenerate Dual Extremum Problem 2.

⁢κ 0 , n ⁡ ( C ) = inf ⁢ E ⁢ { ℓ 0 , n ⁡ ( A n , Y n ) } ⁢ ⁢ if ⁢ ⁢ 𝒫 [ 0 , n ] = 𝒫 [ 0 , n ] D = κ 0 , n ⁡ ( 0 ) .

The last equality follows from fact that randomized control strategies do not incur better performance.

The application example in [0037] illustrates the above statements.

Characterization of Control-Coding Capacity of Dynamical Systems

G. Consider an example control system with control system distribution, and cost function defined by Q _(i)(dy _(i) |y ^(i-1) ,a _(i-L) ^(i)),γ_(i)(a _(i-N) ^(i) ,y ^(i-1)),i=0, . . . ,n  (.1)

where a_(i-L) ^(i)

(a_(i-L), a_(i-L+1), . . . , a^(i)), y^(i-1)=(y⁻¹, y₀, . . . , y_(i-1)), {L, N} are a finite non-negative integers, with an average constraint defined by

${\mathcal{P}_{\lbrack{0,n}\rbrack}^{N}(\kappa)} = {\left\{ {{P_{i}\left( {\left. {da}_{i} \middle| a^{i - 1} \right.,y^{i - 1}} \right)},{i = 0},\ldots\mspace{14mu},{n:{{\frac{1}{n + 1}{E\left( {\sum\limits_{i = 0}^{n}{\gamma_{i}\left( {A_{i - N}^{i},Y^{i - 1}} \right)}} \right)}} \leq \kappa}}} \right\}.}$

Then the characterization of FTH information CC-Capacity is given by

⁢J A n → Y n I ⁡ ( π I , * , κ ) ⁢ sup ⁢ ∑ i = 0 n ⁢ E ⁢ { log ⁡ ( dQ i ⁡ ( · | Y i - 1 , A i - L i ) d ⁢ ⁢ Π i π I ⁡ ( · | Y i - 1 ) ⁢ ( Y i ) ) } , ⁢ ⁢ I = max ⁢ { L , N } ( .2 )

where the distributions are given by

∏ i π I ⁢ ( dy i | y i - 1 ) = ∫ ⁢ Q i ⁡ ( dy i | y i - 1 , a i - L i ) ⊗ π i I ⁡ ( da i | a i - I i - 1 , y i - 1 ) ⊗ P π I ⁡ ( da i - I i - 1 | y i - 1 ) , ( .3 ) P μ π I ⁡ ( da i , dy i ) = μ ⁡ ( dy - 1 ) ⁢ ⊗ i = 0 n ⁢ ( Q i ⁡ ( dy i | y i - 1 , a i - L i ) ⊗ π i I ⁡ ( da i | a i - I i - 1 , y i - 1 ) ) , ⁢ ⁢ i = 0 , … ⁢ , n . ⁢ 𝒫 _ [ 0 , n ] I ⁡ ( κ ) ⁢ = Δ ⁢ { π i 1 ⁡ ( da i | a i - I i - 1 , y i - 1 ) , i = 0 , … ⁢ , n : 1 n + 1 ⁢ E ⁡ ( ∑ i = 0 n ⁢ γ i ⁡ ( A i - N i , Y i - 1 ) ) ≤ κ } ( .4 )

H. Consider an example control system with control system distribution, and cost function defined by Q _(i)(dy _(i) |y _(i-M) ^(i-1) ,a _(i-L) ^(i)),γ_(i)(a _(i-N) ^(i) ,y _(i-K) ^(i-1)),i=0, . . . ,n  (.5)

where {M, K} are finite non-negative integers and the convention is

     y_(i = M)^(i = 1)|_(M = 0) = {Ø} for any i. Then the characterization of FTH information CC-Capacity is given by

⁢J A * → Y * I , M ⁡ ( π I , * , κ ) = sup ⁢ ∑ i = 0 n ⁢ E ⁢ { log ⁡ ( dQ i ⁡ ( · | Y i - M i - 1 , A i - L i ) d ⁢ ∏ i π L ⁢ ( · | Y i - 1 ) ⁢ ( Y i ) ) } ⁢( .6 )

where

_([0,n]) ^(I,K)(κ)

{π_(i) ^(I)(da _(i) |a _(i-1) ^(i-1) ,y ^(i-1)),i=0, . . . ,n:1/n+1E(Σ_(i=0) ^(n)γ_(i)(A _(i-N) ^(i) ,Y _(i-K) ^(i-1)))≤κ}, where I

max{L, N}, and

∏ i π I ⁢ ( dy i | y i - 1 ) = ∫ ⁢ Q i ⁡ ( dy i | y i - M i - 1 , a i - L i ) ⊗ π i 1 ⁡ ( da i | a i - I i - 1 , y i - 1 ) ⊗ P π I ⁡ ( da i - I i - 1 | y i - 1 ) ( .7 ) P μ π I ⁡ ( da i , dy i ) = μ ⁡ ( dy - 1 ) ⁢ ⊗ j = 0 i ⁢ ( Q j ⁡ ( dy j | y j - M i - 1 , a j - L j ) ⊗ π j I ⁡ ( da j | a j - I j - 1 , y j - 1 ) ) , i = 0 , … ⁢ , n . ( .8 )

I. Consider an example control system with control system distribution and cost function defined by Q _(i)(dy _(i) |y _(i-M) ^(i-1) ,a _(i)),γ_(i)(a _(i) ,y _(i-K) ^(i-1)),i=0, . . . ,n.  (.9)

Then the characterization of FTH information CC-Capacity is given by

⁢J A n → Y n J ⁡ ( π 0 , J , κ ) = sup ⁢ ∑ i = 0 n ⁢ E ⁢ { log ⁡ ( dQ i ⁡ ( · | Y i - M i - 1 , A i ) d ⁢ ∏ i π 0 , j ⁢ ( · Y i - J i - 1 ) ⁢ ( B i ) ) } , ⁢⁢J ⁢ = Δ ⁢ max ⁢ { M , K } ( .11 )

where

𝒫 ° [ 0 , n ] J ⁡ ( κ ) ⁢ = Δ ⁢ { π 0 , J ⁡ ( da i | y i - J i - 1 ) , i = 0 , … ⁢ , n : 1 n + 1 ⁢ E ⁡ ( ∑ i = 0 n ⁢ γ i ⁡ ( A i , Y i - K i - 1 ) ) ≤ κ } . ⁢ ⁢ ∏ i π 0 , J ⁢ ( dy i | y i - J i - 1 ) = ∫ Q i ⁡ ( dy i | y i - M i - 1 , a i ) ⊗ π j 0 , J ⁡ ( da i | y i - J i - 1 ) , ⁢ P π π 0 , J ⁡ ( da i , dy i ) = μ ⁡ ( dy - 1 ) ⁢ ⊗ j = 0 i ⁢ ( Q j ⁡ ( dy j | y j - M j - 1 , a j ) ⊗ π j 0 , J ⁡ ( da j | y i - J i - 1 ) ) . ( .11 )

The above characterization means the joint process {(A_(i), Y_(i)): i=0, . . . , n} and the output process {Y_(i): i=0, . . . , n} are J-order Markov process.

J. The characterizations of FTH information CC-Capacity are general and hence, they hold for arbitrary control models, noise distributions, and cost functions.

CC-Capacity of Gaussian-G-RCM-1 and Randomized Strategies

K. Consider an example G-RCM of E. called Gaussian-G-RCM-1 with quadratic cost function defined as follows. Y _(i) =C ^(i-1) Y ^(i-1) +D _(i,i) A _(i) +D _(i,i-1) A _(i-1) +V _(i), Y ⁻¹ =y ⁻¹ ,A ₋₁ =a ₋₁ ,i=0, . . . ,n, P _(V) _(i) _(|V) _(i-1) _(,A) _(i) _(,Y) ₋₁ =P _(V) _(i) ,V _(i) ˜N(0,K _(V) _(i) ),K _(V) _(i) >0, (Y ⁻¹ ,A ₋₁)˜N(0,K _(Y-1) ,A ₋₁),K _(Y-1) ,A ₋₁>0, γ_(i)(a _(i) ,y _(i-1))

a _(i) ,R _(i) a _(i)

+

y _(i-1) ,Q _(i,i-1) y _(i-1)

, (D _(i,i) ,D _(i,i-1))∈

^(p×q)×

^(p×q), R _(i)∈

₊₊ ^(q×q) ,Q _(i,i-1)∈

₊ ^(p×p) ,i=0,β,n.

where

₊ ^(q×q) denotes the set of positive semidefinite q×q matrices, and

₊₊ ^(q×q) their restriction to positive definite, and

⋅,⋅

denotes inner product. The Gaussian-RCM-1 is a Multiple Input Multiple Output (MIMO) control system with memory on past inputs and outputs, is an Infinite Impulse Response (IIR) model, and the cost function is quadratic.

Next, we prepare to compute

(i) the optimal randomized strategy, and

(ii) the FTH CC-Capacity.

From G. the optimal randomized strategy is of the form {π_(i) ^(L)(da_(i)|a_(i-L) ^(i-1),y^(i-1))≡π_(i) ^(l)(da_(i)|a^(i-1),y^(i-1)): i=0, . . . , n}, i.e., L=1. The directed information pay-off is expressed as follows.

$\mspace{79mu}{{{I\left( A^{n}\rightarrow Y^{n} \right)} = {\sum\limits_{i = 0}^{n}\left\{ {{H\left( Y_{i} \middle| Y^{i - 1} \right)} = {H\left( {\left. Y_{i} \middle| A_{i - 1}^{i} \right.,Y^{i - 1}} \right)}} \right\}}},\mspace{79mu}{{H\left( {\left. Y_{i} \middle| Y^{i - 1} \right.,A_{i - 1}^{i}} \right)} = {{H\left( V_{i} \right)} = {\frac{1}{2}{\log\left( {2\pi\; e} \right)}^{p}{{K_{V_{i}}}.}}}}}$

Let {(A_(i) ^(g), Y_(i) ^(g), Z_(i) ^(g)): i=0, . . . , n} denote a jointly Gaussian process. By the maximum entropy property of Gaussian distributions it follows that

${\sum\limits_{i = 0}^{n}{H\left( Y_{i} \middle| Y^{i = 1} \right)}} = {{H\left( Y^{n} \right)} \leq {H\left( Y^{g,n} \right)}}$

and the upper bound is achieved if {(A_(i), Y_(i), Z_(i))=(A_(i) ^(g), Y_(i) ^(g), Z_(i) ^(g)): z=0, . . . , n} and the average constraint is satisfied. Hence, the upper bound is achieved, if the optimal strategies are linear, given as follows.

Randomized Strategy

$\begin{matrix} {{A_{i}^{g} = {U_{i}^{g} + {A_{i,{i - 1}}A_{i - 1}^{g}} + Z_{i}^{g}}},{U_{i}^{g}\overset{\Delta}{=}{\Gamma^{i - 1}Y^{g,{i - 1}}}},} \\ {{\equiv {{g_{i}^{1}\left( Y^{g,{i - 1}} \right)} + {A_{i,{i - 1}}A_{i - 1}^{g}} + Z_{i}^{g}}},} \end{matrix}$ g_(i)¹(y^(i − 1)) = Γ^(i − 1)y^(i − 1), i = 0, …  , n,

-   -   Z_(i) ^(g) is independent of (A^(g,i-1), Y^(g,i-1)), Z^(g,i) is         independent of V^(i), i=0, . . . , n,     -   {Z_(i) ^(g)˜N(0, K_(Z) _(i) ): i=0, 1, . . . , n} is an         independent Gaussian process

for some deterministic matrices {(Γ^(i-1), Λ_(i,i-1)): i=0, . . . , n} of appropriate dimensions. Next, we prepare to compute the directed information pay-off. To this end, we need to compute the conditional entropy H(Y_(i) ^(g)|Y^(g,i-1)), i=0, . . . , n, which means we need to determine the conditional density of Y_(i) ^(g) given Y^(g,i-1) for i=0, . . . , n, using the stochastic control system and strategy. Since the conditional density is characterized by the conditional mean and covariance, we define the quantities Ŷ _(i|i-1)

E^(e) ¹ {Y _(i) ^(g) |Y ^(g,i-1) },Â _(i|i) ΔE ^(e) ¹ {A _(i) ^(g) |Y ^(g,i)}, K _(Y) _(i) _(|Y) _(i-1)

E^(e) ¹ {(Y _(i) ^(g) −Ŷ _(i|i-1))(Y _(i) ^(g) −Ŷ _(i|i-1))^(T) |Y ^(g,i-1)} P _(i|i) =E ^(e) ¹ (A _(i) ^(g) −Â _(i|i))(A _(i) ^(g) −Â _(i|i))^(T) ,i=0, . . . ,n.

From above and using the independent properties of the noise process then Â _(i|i)=Λ_(i,i-1) Â _(i-1|i-1+1) U _(i) ^(g)+Δ_(i|i-1)(Y _(i) ^(g) −Ŷ _(i|i-1)), Ŷ _(i|i-1) =C ^(i-1) Y ^(g,i-1) +D _(i,i) Y _(i) ^(g)+Λ _(i,i-1) Â _(i-1|i-1), K _(Y) _(i) _(|Y) _(i-1) =Λ _(i,i-1) P _(i-1|i-1) Λ _(i,i-1) ^(T) +D _(i,i) K _(Z) _(i) D _(i,i) ^(T) +K _(V) _(i) ,i=0, . . . ,n

where (Â_(−1|−1), P_(−1|−1)) are initial data and Λ _(i,i-1)

D_(i,i)Λ_(i,i-1) +D _(i,i-1) ,i=0, . . . ,n, P _(i|i)=Λ_(i,i-1) P _(i-1|i-1)Λ_(i,i-1) ^(T) +K _(Z) _(i) −(K _(Z) _(i) D _(i,i) ^(T)+Λ_(i,i-1) P _(i-1|i-1) Λ _(i,i-1) ^(T))·Φ_(i|i-1)(K _(Z) _(i) D _(i,i) ^(T)+Λ_(i,i-1) P _(i-1|i-1) Λ _(i,i-1) ^(T))^(T), Φ_(i|i-1)

[D _(i,i) K _(Z) _(i) D _(i,i) ^(T) +K _(V) _(i) =Λ _(i,i-1) P _(i-1|i-1) Λ _(i,i-1) ^(T)]⁻¹, Δ_(i|i-1)

(K _(Z) _(i) D _(i,i) ^(T)+Λ_(i,i-1) P _(i-1|i-1) Λ _(i,i-1) ^(T))Φ_(i|i-1)

The innovations process denoted by {v^(e) ¹ : i=0, . . . , n} is an orthogonal process, independent of {g_(i) ¹(⋅): i=0, . . . , n}, and satisfies the following identities.

${v_{i}^{e^{1}}\overset{\Delta}{=}{{Y_{i}^{g} - {\hat{Y}}_{{i|i} = 1}} = {{{{\overset{\_}{A}}_{i,{i = 1}}\left( {A_{i = 1}^{g} = {\hat{A}}_{{i - 1} \mid {i - 1}}} \right)} + {D_{i,i}Z_{i}^{g}} + V_{i}} = \left. v_{i}^{e^{2}} \middle| {}_{g^{i} = 0}{\equiv v_{i}^{0}} \right.}}},{\left. v_{i}^{0} \right.\sim{N\left( {0,K_{Y_{i}|Y^{i - s}}} \right)}},{i = 0},\ldots\mspace{14mu},n$

where {v_(i) ⁰: i=0, . . . , n} indicates that the innovations process is independent of the strategy {g_(i) ¹(⋅): i=0, . . . , n}. From the above equations, since the conditional covariance K_(Y) _(i) _(|Y) _(i-1) is independent of Y^(g,i-1) for i=0, . . . , n, the conditional distribution is P^(e) ¹ (Y_(i) ^(g)≤y_(i)|Y^(g,i-1))˜N(Ŷ_(i|i-1), K_(Y) _(i) _(|Y) _(i-1) ), i=0, . . . , n. Applying the above two observations we obtain

$\mspace{79mu}{{I\left( A^{g,n}\rightarrow Y^{g,n} \right)} = {{\frac{1}{2}{\sum\limits_{i = 0}^{n}{\log\frac{K_{Y_{i}|Y^{i - 1}}}{K_{V_{i}}}}}} \equiv {\sum\limits_{i = 0}^{n}\left\{ {{H\left( v_{i}^{0} \right)} - {H\left( V_{i} \right)}} \right\}}}}$

Next, we give the the closed form expressions of the optimal randomized control strategies, and the FTH CC-Capacuty.

CC-Capacity and Randomized Strategy

Consider the Gaussian-RCM-1. The following hold.

(a) FTH Information CC-Capacity. The joint process {(A_(i), Y_(i))=(A_(i) ^(g), Y_(i) ^(g)), i=0 . . . , n}, is jointly Gaussian and satisfies the following equations.

$\begin{matrix} {{A_{i}^{g} = {{e_{i}^{1}\left( {Y^{s,{i - 1}},A_{i - 1}^{g}} \right)} + Z_{i}^{g}}},{i = 0},\ldots\mspace{14mu},n,} \\ {{= {U_{i}^{g} + {A_{i,{i - 1}}A_{i - 1}^{g}} + Z_{i}^{g}}},} \end{matrix}$ ${U_{i}^{g} = {{g_{i}^{1}\left( Y^{g,{i - 1}} \right)} = {\Gamma^{i - 1}Y^{g,{i - 1}}}}},{Y_{i}^{g} = {{C^{i - 1}Y^{g,{i - 1}}} + {{\overset{\_}{A}}_{i,{i - 1}}A_{i - 1}^{g}} + {D_{i,i}U_{i}^{g}} + {D_{i,i}Z_{i}^{g}} + V_{i}}},$

-   -   i) Z_(i) ^(g) is independent of (A^(g,i-1), Y^(g,i-1)), i=0, . .         . , n,     -   ii) Z^(g,i) is independent of V^(i), i=0, . . . , n,     -   iii) {Z_(i) ^(g)˜N(0, K_(Z) _(i) ): i=0, 1, . . . , n} is an         independent Gaussian process,

E^(e¹){γ_(i)(A_(i)^(g), Y_(i − 1)^(g))} = E^(e¹){⟨U_(i)^(g), R_(i)U_(i)^(g)⟩ + 2⟨A_(i, i − 1)Â_(i − 1|i − 1), R_(i)U_(i)^(g)⟩ + ⟨A_(i, i − 1)Â_(i − 1|i − 1), R_(i)A_(i, i − 1)Â_(i, i − 1|i − 1)⟩ + tr(K_(Z_(i))R_(i)) + tr(A_(i, i − 1)^(T)R_(i)A_(i, i − 1)P_(i − 1|i − 1)) + ⟨Y_(i − 1)^(g), Q_(i)Y_(i − 1)^(s)⟩}.

The FTH CC-Capacity given by

⁢J A n → Y n 1 ⁡ ( e 1 , * ) = sup ⁢ 1 2 ⁢ ∑ i = 0 n ⁢ log ⁢  K Y i | Y i - 1   K V i 

and the average constraint set is defined by

${{\overset{\_}{\mathcal{P}}}_{\lbrack{0,n}\rbrack}^{1}(\kappa)}\overset{\Delta}{=}{\left\{ {{{e_{i}^{1}( \cdot )}\overset{\Delta}{=}\left( {{g_{i}^{1}\left( {\cdot {, \cdot}} \right)},A_{i,{i - 1}},K_{Z_{i}}} \right)},{i = 0},\ldots\mspace{14mu},{n:{{\frac{1}{n + 1}{\sum\limits_{i = 0}^{n}{E^{e^{2}}\left( {\gamma_{i}\left( {A_{i}^{g},Y^{g,{i - 1}}} \right)} \right)}}} \leq \kappa}}} \right\}.}$

(b) Decentralized Separation of Randomized Strategy into Controller and Encoder: The optimal strategy denoted by {e^(1,*)(⋅)≡(g_(i) ^(1,*)(⋅), Λ_(i,i-1)*, K_(Z) _(i) *): i=0, . . . , n} is the solution of the dual optimization problem

${\kappa_{0,n}(C)}\overset{\Delta}{=}{\inf\limits_{{({{g_{i}^{1}{( \cdot )}},A_{i,{i - 1}},K_{Z_{i}}})},{i = 0},\ldots\mspace{14mu},{n:{{\frac{1}{2}{\sum\limits_{i = 0}^{n}{\log\frac{K_{Y \mid Y^{i - 1}}}{K_{V_{i}}}}}} \geq {{({n + 1})}C}}}}E^{e^{2}}{\left\{ {\sum\limits_{i = 0}^{n}{\gamma_{i}\left( {A_{i}^{g},Y_{i - 1}^{g}} \right)}} \right\}.}}$

Moreover, the following decentralized separation holds.

(i) The optimal strategy {g_(i) ^(1.*)(⋅): i=0, . . . , n} is the solution of the optimization problem

$\inf\limits_{{{g_{i}^{1}{( \cdot )}};{i = 0}},\ldots\mspace{14mu},n}E^{e^{2}}\left\{ {\sum\limits_{i = 0}^{n}{\gamma_{i}\left( {A_{i}^{g},Y_{i - 1}^{g}} \right)}} \right\}$

for a fixed {Λ_(i,i-1), K_(Z) _(i) : i=0, . . . , n}.

(ii) The optimal strategy {Λ_(i,i-1)*, K_(Z) _(i) *: i=0, . . . , n} is the solution of κ_(0,n)(C) for {g_(i) ¹(⋅)=g_(i) ^(1,*)(⋅): i=0, . . . , n} determined from (i).

(c) Optimal Strategies of Controller and Encoder. Suppose in (a), Y_(i) ^(g) is replaced by Y _(i) ^(g) =C _(i,i-1) Y _(i-1) ^(g)+Λ _(i,i-1) A _(i-1) ^(g) +D _(i,i) U _(i) ^(g) +D _(i,i) Z _(i) ^(g) +V _(i) ,i=0, . . . ,n.

Any candidate of the control strategy {g_(i) ¹(Y_(g,i-1)): i=0, . . . , n} is of the form

$\begin{matrix} {{{g_{i}^{1}\left( Y^{g,{i - 1}} \right)}\overset{\Delta}{=}{{\Gamma_{i,{i - 1}}^{1}Y_{i - 1}^{g}} + {\Gamma_{i,{i - 1}}^{2}{\hat{A}}_{{i - 1}|{i - 1}}}}},} \\ {{\equiv {{\overset{\_}{\Gamma}}_{i,{i - 1}}{\overset{\_}{Y}}_{i - 1}^{g}}},{{\overset{\_}{Y}}_{i - 1}^{g} = \begin{bmatrix} Y_{i - 1}^{g} \\ {\hat{A}}_{{i - 1}|{i - 1}} \end{bmatrix}},{i = 0},\ldots\mspace{14mu},{n.}} \end{matrix}$

Define the augmented system

${{\overset{\_}{Y}}_{i}^{g} = {{{\overset{\_}{F}}_{i,{i - 1}}{\overset{\_}{Y}}_{i - 1}^{g}} + {{\overset{\_}{B}}_{i,{i - 1}}U_{i}^{g}} + {{\overset{\_}{G}}_{i,{i - 1}}v_{i}^{e^{2}}}}},{{\overset{\_}{F}}_{i,{i - 1}}\overset{\Delta}{=}\begin{bmatrix} C_{i,{i - 1}} & {\overset{\_}{A}}_{i,{i - 1}} \\ 0 & A_{i,{i - 1}} \end{bmatrix}},{{\overset{\_}{B}}_{i,{i - 1}}\overset{\Delta}{=}\begin{bmatrix} D_{i,i} \\ I \end{bmatrix}},{{\overset{\_}{G}}_{i,{i - 1}}\overset{\Delta}{=}\begin{bmatrix} I \\ \Delta_{i|{i - 1}} \end{bmatrix}},{i = 0},\ldots\mspace{14mu},n$

and average cost

E e 1 ⁢ { ∑ i = 0 n ⁢ γ i ⁢ ( A i g , Y i - 1 g ) } ≡ ⁢ E ⁢ e 1 ⁢ { ∑ i = 0 n ⁢ γ _ i ⁡ ( U i g , Y i - 1 g ) } = Δ ⁢ ⁢ E e 1 ⁢ { ∑ i = 0 n ⁢ ( [ Y _ i - 1 g U i g ] T ⁡ [ M _ i , i - 1 L _ i , i - 1 L _ i , i - 1 N _ i , i - 1 ] ⁡ [ Y _ i - 1 g U i g ] + ⁢ tr ⁡ ( K Z i ⁢ R i ) + tr ⁡ ( A i , i - 1 T ⁢ R i ⁢ A i , i - 1 ⁢ P i - 1 | i - 1 ) ) } . $\mspace{79mu}{{{\overset{\_}{M}}_{i,{i - 1}}\overset{\Delta}{=}\begin{bmatrix} Q_{i,{i - 1}} & 0 \\ 0 & {A_{i,{i - 1}}^{T}R_{i}A_{i,{i - 1}}} \end{bmatrix}},\mspace{79mu}{{\overset{\_}{L}}_{i,{i - 1}}\overset{\Delta}{=}\begin{bmatrix} 0 \\ {A_{i,{i - 1}}^{T}R_{i}} \end{bmatrix}},{{\overset{\_}{N}}_{i,{i - 1}}\overset{\Delta}{=}{R_{i}.}}}$

Then the following hold.

(1) For a fixed {Λ_(i,i-1), K_(Z) _(i) : i=0, . . . , n} the optimal strategy {U_(i) ^(g,*)=g_(i) ^(1,*)(Y^(g,i-1)): i=0, . . . , n} is the solution of the partially observable classical stochastic optimal control problem

⁢J 0 , n ⁡ ( g 1 , * ⁡ ( · ) , A , K Z ) = inf ⁢ E e 1 ⁢ { ∑ i = 0 n ⁢ γ _ i ⁡ ( U i g , Y _ i = 1 g ) }

where {Y _(i) ^(g): i=0, . . . , n} satisfy the above recursion. Moreover, the optimal strategy {U_(i) ^(g,*)=g_(i) ^(1,*)(Y^(g,i-1))≡g_(i) ^(1,*)(Y ^(g,i-1)): i=0, . . . , n} is given by the following equations.

$\mspace{20mu}{{{g_{i}^{1,*}\left( {\overset{\_}{y}}^{i - 1} \right)} = {{\overset{\_}{\Gamma}}_{i,{i = 1}}{\overset{\_}{y}}^{i = 1}}},{i = 0},\ldots\mspace{11mu},{n - 1},{= {{- \left( {{\overset{\_}{N}}_{i,{i - 1}} + {{\overset{\_}{B}}_{i,{i - 1}}{\sum{\left( {i + 1} \right){\overset{\_}{B}}_{i,{i - 1}}}}}} \right)^{- 1}}{\overset{\_}{B}}_{i,{i - 1}}^{T}{\sum{\left( {i + 1} \right){\overset{\_}{F}}_{i,{i - 1}}{\overset{\_}{y}}^{i - 1}}}}},\mspace{20mu}{{g_{n}^{1,*}\left( {\overset{\_}{y}}^{n - 1} \right)} = 0},}$

where the symmetric positive semidefinite matrix {Σ(i): i=0, . . . , n} satisfies a matrix difference Riccati equation, for i=0, . . . , n−1, equation Σ(i)= F _(i,i-1) ^(T)Σ(i+1) F _(i,i-1)−( F _(i,i-1) ^(T)Σ(i+1) B _(i,i-) +L _(i,i-1))·( N _(i,i-1) +B _(i,i-1) ^(T)Σ(i+1) B _(i,i-1))⁻¹( B _(i,i-1) ^(T)Σ_(i,i-1) F _(i,i-1) +L _(i,i-1) ^(T))^(T) +M _(i,i-1),Σ(n)= M _(n,n-1) and the optimal pay-off is given by

${J_{0,n}\left( {{g^{1,*}( \cdot )},A,K_{Z}} \right)} = {\sum\limits_{j = 0}^{n}{\left\{ {{{{tr}\left( {\Phi_{j \mid {j - 1}}{\sum(j)}} \right)} + {{tr}\left( {{K_{Z_{j}}R_{j}} + {{tr}\left( {A_{j,{j - 1}}^{T}R_{j}A_{j,{j - 1}}P_{{j - 1}|{j - 1}}} \right)} + {{tr}\left( {\Delta_{j|{j - 1}}K_{Y_{j}|Y_{j - 1}}\Delta_{{jj} - 1}{\sum(j)}} \right)}} \right\}} + {\overset{\_}{Y}}_{{- 1}|{- 1}}},{\sum{(0){\overset{\_}{Y}}_{{- 1}|{- 1}}}}} \right).}}$

(2) The optimal strategies {(Λ_(i,i-1)*, K_(Z) _(i) *): i=0, . . . , n} are the solutions of the optimization problem

$\mspace{79mu}{{\kappa_{0,n}(C)} = {\inf\limits_{{({A_{i,{1 - ɛ}},K_{x_{i}}})},{A = 0},{{\ldots\mspace{14mu} n};{{\frac{1}{2}{\sum\limits_{i = 0}^{n}{\log\frac{K_{Y_{i} \mid {{Yi} - 1}}}{K_{V_{i}}}}}} \geq {{({n + 1})}C}}}}{\left\{ {J_{0,n}\left( {{g^{1,*}( \cdot )},A,K_{Z}} \right)} \right\}.}}}$

The above is a decentralized separation principle, and (1) and (2) are Person-by-Person Optimality statements of {g_(i) ¹(⋅): i=0, . . . ,} and {Λ_(i,i-1), K_(Z) _(i) : i=0, . . . , n}.

CC-Capacity of Gaussian-G-RCM-2 and Randomized Strategies

L. Consider an example G-RCM of H. called Gaussian-G-RCM-2 with quadratic cost function defined as follows. Y _(i) =C _(i,i-1) Y _(i-1) +D _(i) A _(i) +V _(i) ,Y ₋₁ =y ₋₁, P _(V) _(i) _(|V) _(i-1) _(,A) _(i) _(,Y) ₋₁ =P _(V) _(i) ,V _(i) ˜N(0,K _(V) _(i) ),K _(V) _(i) >0, γ_(i)(a _(i) ,y _(i-1))

a _(i) ,R _(i) a _(i)

+

y _(i-1) ,Q _(i,i-1) y _(i-1)

.

By I. the characterization if FTH-DI information CC-Capacity is

J A n → Y n ⁡ ( π * , κ ) = ⁢ sup ° [ 0 , n ] ⁡ ( κ ) ⁢ ∑ i = 0 n ⁢ ∫ log ⁡ ( dQ i ⁡ ( · | y i - 1 , a i ) d ⁢ ∏ i * ⁢ ( · | y i - 1 ) ⁢ ( y i ) ) ⁢ P * ⁡ ( dy i , y i = 1 , da i ) ≡ ⁢ sup ⁢ ∑ i = 0 n ⁢ I ⁡ ( A i ; Y i | Y i - 1 ) ⁢ ∏ i * ⁢ ( dy i | y i - 1 ) = ∫ λ i ⁢ Q i ⁡ ( da i | y i - 1 , a i ) ⊗ π i ⁡ ( da i | y i - 1 ) 𝒫 Δ [ 0 , n ] ⁡ ( κ ) ⁢ = Δ ⁢ { π i ⁡ ( da i | y i - 1 ) , i = 0 , … ⁢ , n : 1 n + 1 ⁢ ∑ i = 0 n ⁢ E n ⁢ { 〈 A i , R i ⁢ A i 〉 + 〈 Y i - 1 , Q i , i - 1 ⁢ Y i - 1 〉 } ≤ κ } .

The optimal randomized strategy of above characterization of FTH information CC-Capacity is now computed, using several steps.

(a) Gaussian Properties of Characterization of FTH Information CC-Capacity. Using dynamic programming or the maximum entropy property of processes with fixed second moments, the optimal strategies are Gaussian denoted by {π_(i) ^(g)(da_(i)|y_(i-1)): i=0, . . . , n}∈P̊_([0,n])(κ), and the joint process is jointly Gaussian denoted by {(A_(i), Y_(i))≡(A_(i) ^(g), Y_(i) ^(g)): i=0, . . . , n}.

(b) Realization of Optimal Strategies. Since {(A_(i) ^(g), Y_(i) ^(g)): i=0, . . . , n} is jointly Gaussian, strategies from the set P̊_([0,n])(κ), can be realized by linear and Gaussian randomized strategies defined by the set

${\mathcal{P}_{\lbrack{0,n}\rbrack}(\kappa)} = \left\{ {{A_{i}^{g} = {e_{i}\left( {Y_{i = 1}^{g},Z_{i}^{g}} \right)}},{{e_{i}\left( {Y_{i = 1}^{g},Z_{i}^{g}} \right)} = {{{g_{i}\left( Y_{i = 1}^{g} \right)} + Z_{i}^{g}} = {{\Gamma_{i,{i - 1}}Y_{i = 1}^{g}} + Z_{i}^{g}}}},{Z_{i}^{g}\bot Y^{g,{i = 1}}},{\left\{ {{{Z_{i}^{g}:i} = 0},\ldots\mspace{14mu},n} \right\}\mspace{14mu}{independent}\mspace{14mu}{process}},{\left. Z_{i}^{g} \right.\sim{N\left( {0,K_{Z_{i}}} \right)}},{K_{Z_{i}} \in S_{+}^{q \times q}},{i = 0},\ldots\mspace{14mu},{n:{{\frac{1}{n + 1}E^{e}\left\{ {\sum\limits_{i = 0}^{n}\left\lbrack {\left\langle {A_{i}^{g},{R_{i}A_{i}^{g}}} \right\rangle + \left\langle {Y_{i = 1}^{g},{Q_{i,{i = 1}}Y_{i = 1}^{g}}} \right\rangle} \right\rbrack} \right\}} \leq \kappa}}} \right\}$

where ⋅⊥⋅ means the processes are independent.

Moreover, the characterization of the FTH information CC-Capacity is

⁢J A n → Y n ⁡ ( π 3 , * , κ ) = sup ⁢ { 1 2 ⁢ ∑ i = 0 n ⁢ ln ⁢ | D i ⁢ K Z i ⁢ D i T + K V i | | K V i | } . (c) Dual Role of Randomized Control Strategies. Since the optimal control strategies admit the decomposition A _(i) ^(g)=Γ_(i,i-1) Y _(i-1) ^(g) +Z _(i) ^(g) ≡g _(i)(Y _(i-1) ^(g))+Z _(i) ^(g) ,i=0, . . . ,n

then we have the following: (i) the feedback control law or strategy {g_(i)≡Γ_(i,i-1): i=0, . . . , n} is responsible for controlling the output process {Y_(i) ^(g): i=0, . . . , n}, and (ii) the orthogonal innovations process {Z_(i) ^(g): i=0, . . . , n} is responsible for communicating new information to the output process, both chosen to maximize J_(A) _(n) _(→Y) _(n) (π^(g,*),κ).

-   -   (d) Decentralize Separation Principle of Optimal Randomized         Strategies. Let {(A_(i) ^(g,*), T_(i) ^(g,*), Z_(i) ^(g,*)):         i=0, . . . , n} denote the optimal joint process corresponding         to (c) above. The cost-to-go C_(i):         _(i-1)         (corresponding to J_(A) _(n) _(→Y) _(n) (π^(g,*), κ)), from time         “i” to terminal time “n” satisfies the dynamic programming         recursions

C n ⁡ ( y n - 1 ) = sup ⁢ { 1 2 ⁢ log ⁢  D n ⁢ K Z n ⁢ D n T + K V n   K V n  - str ⁡ ( R n ⁢ K Z n ) - s ⁡ [ 〈 u b , R n ⁢ u n 〉 + 〈 y n - 1 , Q n , n - 1 ⁢ y n - 1 〉 ] + s ⁡ ( n + 1 ) ⁢ κ } , ⁢ C i ⁡ ( y i - 1 ) = sup ⁢ { 1 2 ⁢ log ⁢  D i ⁢ K Z i ⁢ D i T + K V i   K V i  - str ⁡ ( R i ⁢ K Z i ) - s ⁡ [ 〈 i i , R i ⁢ u i 〉 + 〈 y i - 1 , Q i , i - 1 ⁢ y i - 1 〉 ] + E ⁢ { C i + 1 ⁡ ( Y i g , * ) | Y i - 1 g , * = y i - 1 } } , i = 0 , … ⁢ , n - 1

where s≥0 is the Lagrange multiplier associated with the average constraint. The solution of the dynamic programming equations is given by the following equations. C _(i)(y _(i-1))={−s(y _(i-1) ,P(i)y _(i-1))+r(i)},i=0, . . . ,n

where {r(i): i=0, . . . , n−1} satisfies the recursions

${r(i)} = {{r\left( {i + 1} \right)}\sup\limits_{{\epsilon S}_{+}^{q \times q}}\left\{ {{\frac{1}{2}\log\frac{{{D_{i}K_{Z_{i}}D_{i}^{T}} + K_{V_{i}}}}{K_{V_{i}}}} - {{str}\left( {{P\left( {i + 1} \right)}\left\lbrack {{D_{i,i}K_{Z_{i}}D_{i}^{T}} + K_{V_{i}}} \right\rbrack} \right)} - {{str}\left( {R_{i}K_{Z_{i}}} \right)}} \right\}}$ ${r(n)} = {\sup\limits_{{\epsilon S}_{+}^{q \times q}}\left\{ {{\frac{1}{2}\log\frac{{{D_{n}K_{Z_{n}}D_{n}^{T}} + K_{V_{n}}}}{K_{V_{n}}}} - {{str}\left( {R_{n}K_{Z_{n}}} \right)} + {{s\left( {n + 1} \right)}\kappa}} \right\}}$

and {P(i): i=0, . . . , n} is a solution of the Riccati difference matrix equation P(i)=C _(i,i-1) ^(T) P(i+1)C _(i,i-1) +Q _(i,i-1) −C _(i,i-1) ^(T) P(i+1)D _(i)(D _(i,i) ^(T) P(i+1)D _(i) +R _(i,i))⁻¹(C _(i,i-1) ^(T) P(i+1)D _(i))^(T) ,i=0, . . . ,n−1, P(n)=Q _(n,n-1).

The optimal randomized control strategy is given by A _(i) ^(g,*) =g _(i)*(Y _(i-1) ^(g,*))+Z _(i) ^(g,*) ,i=0, . . . ,n

where its random part {Z_(i) ^(g,*): i=0, . . . , n} is the solution to the above recursions, and its deterministic part is given by

$\begin{matrix} \begin{matrix} {{g_{i}^{*}\left( y_{t - 1} \right)} = {{- \left( {{D_{i}^{T}{P\left( {i + 1} \right)}D_{i}} + R_{i}} \right)^{- 1}}D_{i}^{T}{P\left( {i + 1} \right)}C_{i,{i - 1}}y_{i - 1}}} \\ {{\equiv {\Gamma_{i,{i - 1}}^{*}y_{i - 1}}},{i = 0},\ldots\mspace{14mu},{n - 1},} \end{matrix} & \; \\ {{g_{n}^{*}\left( y_{n - 1} \right)} = 0.} & \; \end{matrix}$

The corresponding covariance K_(Y) _(i)

E{Y_(i) ^(g)(Y_(i) ^(g))^(T)}, i=0, 1, . . . , n, is given by the following equation. K _(Y) _(i) =(C _(i,i-1) +D _(i)Γ_(i,i-1)*)K _(Y) _(i 1) (C _(i,i-1) +D _(i)Γ_(i,i-1)*)^(T) +D _(i) K _(Z) _(i) D _(i) ^(T) +K _(V) _(i) ,i=0, . . . ,n,K _(Y) ₋₁ =given.

The Lagrange multiplier s≥0 are found from the problem inf _(s≥0) {−s

y ₋₁ ,P(0)y ₋₁

+r(0)}

The characterization of the FTH-DI extremum problem is given by J _(A) _(n) _(→Y) _(n) (π^(g,*),κ)=−s∫ _(Y) ₋₁

y ₋₁ ,P(0)y ₋₁

P _(Y) ₋₁ (dy ₋₁)+r(0).

where s is the value found above.

(e) Water-filling Solution of Encoder Part of Randomized Strategy. The above solution illustrates the decentralized separation between the computation of the optimal deterministic part {g*_(i)(y_(i-1)): i=0, . . . , n} and the optimal random part {K*_(Z) _(i) : i=0, . . . , n} (covariance of innovations process) of the randomized control strategy, the later is found from a sequential water filling problem of recursions r_(i)(⋅) that depends on the solution of the matrix Riccati difference equation. Clearly, optimal randomized control strategies have a dual role; to control the controlled process, precisely as in stochastic optimal control theory, of Gaussian control systems with quadratic pay-off, and to transmit new information via its random part {Z_(i) ^(g): i=0, . . . , n}. This means, information data are mapped into the optimal innovations process with covariance {K*_(Z) _(i) : i=0, . . . , n}, and then transmitted over the stochastic system, which acts as a communication channel.

(f) Connection to Stochastic Optimal Control Problem. From the above solution, we can recover, as a degenerate case, the optimal strategies of Gaussian control problems with quadratic pay-off as follows.

The dual of ( ) is given by

$\begin{matrix} {{\kappa_{0,n}(C)}\overset{\Delta}{=}\inf\limits_{{({\Gamma_{i,{i - 1}},K_{Z_{i}}})},{i = 0},\ldots\mspace{14mu},{n:{{\frac{1}{2{({n + 1})}}{\sum\limits_{i = 0}^{n}{\ln\frac{{{D_{i}K_{Z_{i}}D_{i}^{T}} + K_{V_{i}}}}{K_{V_{i}}}}}} \geq C}}}} \\ {\left\{ {\sum\limits_{i = 0}^{n}{E^{g}\left( {\left( {A_{i}^{g},{R_{i}A_{i}^{g}}} \right) + \left( {Y_{i - 1}^{g},{Q_{i,{i - 1}}Y_{i - 1}^{g}}} \right)} \right)}} \right\}} \\ {= \inf\limits_{{A_{i}^{g} = {g_{i}{(Y_{i - 1}^{g})}}},{i = 0},\ldots\mspace{14mu},n}} \\ {{E^{g}\left\{ {\sum\limits_{i = 0}^{n}\left( {\left( {A_{i}^{g},{R_{i}A_{i}^{g}}} \right) + \left( {Y_{i - 1}^{g},{Q_{i,{i - 1}}Y_{i - 1}^{g}}} \right)} \right)} \right\}},} \\ {{{{if}\mspace{14mu} K_{Z_{i}}} =},{i = 0},\ldots\mspace{14mu},{n.}} \end{matrix}$

The second identity holds if randomized control strategies are restricted to deterministic strategies. Hence, J _(A) _(n) _(→Y) _(n) (#^(g,*),κ)=0 if K* _(Z) _(i) =0,i=0, . . . ,n.

and the degenerate optimization problem is the stochastic optimal control of Gaussian control systems with quadratic pay-off.

CC-Capacity: Per Unit Time of FTH Information CC-Capacity

M. Consider the Gaussian-G-RCM-2 of L. Here, we compute the capacity of the control system, i.e., given by the per unit time limit of the FTH information CC-Capacity. Suppose the control system is time-invariant, with {C_(i,i-1)=C, D_(i)=D, K_(V) _(i) , =K_(V), R_(i)=R, i=0, . . . , n, Q_(i,i-1)=Q, i=0, . . . , n−1, Q_(n,n-1)=M}, and the following conditions hold.

-   -   i) the pair (C, D) is stabilizable,     -   ii) the pair (G, C) is detectable, Q=G^(T)G, G∈         ₊ ^(p×p).

Then the CC-Capacity of the control system is per unit time limit of the characterization of FTH information CC-Capacity, given by

$\begin{matrix} {{C(\kappa)} = {{J_{A^{\infty}\rightarrow Y^{\infty}}\left( {\pi^{g,*},\kappa} \right)}\overset{\Delta}{=}{\lim\limits_{n\rightarrow}{\frac{1}{n + 1}{J_{A^{n}\rightarrow Y^{n}}\left( {\pi^{g,*},\kappa} \right)}}}}} \\ {= {{J^{*}\left( s^{*} \right)} = {\inf\limits_{s \geq 0}{J^{*}(s)}}}} \end{matrix}$      where ${J^{*}(s)} = {\sup\limits_{{\epsilon\mathbb{S}}_{+}^{q \times q}}\left\{ {{\frac{1}{2}\log\frac{{{{DK}_{Z}D^{T}} + K_{V}}}{K_{V}}} + {s\;\kappa} - {{str}\left( {RK}_{Z} \right)} - {{str}\left( {P\left\lbrack {{{DK}_{Z}D^{T}} + K_{V}} \right\rbrack} \right)}} \right\}}$      g^(∞, *)(y) = −(D^(T)PD + R)⁻¹D^(T)PCy ≡ Γ^(∞, *)y.     P = C^(T)PC + Q − C^(T)PD(D^(T)PD + R)⁻¹(C^(T)PD)^(T),     spec(C + D Γ^(∞, *)) ⋐ 𝔻_(o)

where spec(⋅) denotes the set of eigenvalues and

₀

{c∈C:|c|<1} is the unit disc of the set of complex numbers C . The Lagrange multiplier s*(κ) can be found from the average constraint tr(RK _(Z))+tr(P[[DK _(Z) D ^(T) +K _(V)])≤κ.

Thus, the predictable part of the optimal randomized control strategy g_(∞,*)(y) ensures existence of a unique invariant distribution P_(Y) ^(g∞,*)(dy) of the optimal output process {Y*_(i): i=0, . . . , n}corresponding to (g^(∞,*)(⋅), K*_(Z)), i.e., stability of the closed loop system, and hence C(K) is operational, i.e., it is the CC-capacity of the control system.

N. Consider an example Gaussian-G-RCM-2 of L., with parameters p=q=1, R=1, Q=0 and C, D) arbitrary. For these choices of parameters we have the following.

$\begin{matrix} {\left( {\Gamma^{\infty,*},K_{Z}^{*}} \right) = \left\{ {{\begin{matrix} {\left( {0,\kappa} \right),} & {\kappa \in \left\lbrack {0,\infty} \right)} & {{{if}\mspace{14mu}{C}} < 1} \\ {\left( {{- \frac{C^{2} - 1}{CD}},\frac{{D^{2}\kappa} + {K\;{\gamma\left( {1 - C^{2}} \right)}}}{C^{2}D^{2}}} \right),} & {{\kappa \in \left\lbrack {\kappa_{\min},\infty} \right)},{\kappa_{\min} = \frac{\left( {C^{2} - 1} \right)K\;\gamma}{D^{2}}}} & {{{if}\mspace{14mu}{C}} > 1} \\ {\left( {{- \frac{C^{2} - 2}{CD}},0} \right),} & {\kappa \in \left\lbrack {0,\kappa_{\min}} \right\rbrack} & {{{if}\mspace{14mu}{C}} > 1} \end{matrix}\mspace{79mu}{and}\mspace{79mu}{C(\kappa)}} = \left\{ \begin{matrix} {\frac{1}{2}\ln\frac{{D^{2}\kappa} + K_{V}}{K_{V}}} & {{{{if}\mspace{14mu}{C}} < 1},{i.e.},{K_{Z}^{*} = \kappa}} \\ {\frac{1}{2}\ln\frac{{D^{2}K_{Z}^{*}} + K_{V}}{K_{V}}} & {{{{if}\mspace{14mu}{C}} > 1},{\kappa \in \left\lbrack {\kappa_{\min,}\infty} \right)}} \\ 0 & {{{{if}\mspace{14mu}{C}} > 1},{\kappa \in {\left\lbrack {0,\kappa_{\min}} \right\rbrack.}}} \end{matrix} \right.} \right.} & ({.12}) \end{matrix}$

Clearly, if |C|<1, i.e., stable, the deterministic part of the strategy is zero, i.e., Γ^(∞,*)=0. The capacity formulae C(κ), illustrates that there are multiple regimes, depending on whether the control system is stable, that is, |C|<1 or unstable |C|>1. Moreover, for unstable control systems |C|>1, the optimal pay-off is zero, unless the power level r exceeds the critical level κ_(min).

Encoder Design for CC-Capacity of Gaussian-G-RCM-2

O. Consider an example Gaussian-G-RCM-2 of L.

We illustrate via an example of an information process or tracking process, as shown in FIG. 1, how to design the encoder.

Consider a Gaussian information process {X_(i): i=0, 1, . . . , n} taking values in

^(q) (which is to be encoded by the optimal randomized strategy) and described by a Gaussian Linear State Space Model (G-LSSM), as follows. X _(i+1) =A _(i) X _(i) +G _(i) W _(i) ,X ₀ =x,i=0, . . . ,n−1

where {W_(i)˜(0, K_(W) _(i) ): i=0, . . . , n−1} are

=

^(k)—valued zero mean Gaussian processes, independent of the noises driving the Gaussian-G-RCM-2 of the control system, i.e., {(A_(i), V_(i)): i=0, 1, . . . , n}, and the initial state X₀ is either fixed, i.e., X₀=x or its distribution is fixed and Gaussian, i.e., P_(X) ₀ (dx)˜N(0, K_(X) ₀ ), and independent of {W_(i): i=0, . . . , n}. Next, we illustrate how to transform the optimal randomized strategy of L. into an controller-encoder, and a mean-square error decoder, based on a variation of the Kalman filter. Consider the optimal randomized strategy of L., specifically, {(Γ*_(i,i-1), K*_(Z) _(i) ): i=0, 1, . . . , n−1,}, with corresponding optimal distribution {π_(i) ^(g,*)(da_(i)|y_(i-1)): i=0, . . . , n} and joint process {(A_(i) ^(g,*), Y_(i) ^(g,*)): i=0, . . . , n}. Define the filter estimate and conditional covariance, for i=0, . . . , n, by {circumflex over (X)} _(i|i-1)

E{X_(i) |Y ^(g,*,i-1)}, Σ_(i|i-1)

E{(X _(i) −{circumflex over (X)} _(i|i-1))(X _(i) −{circumflex over (X)} _(i|i-1))^(T) |Y ^(g,*,i-1)}.

(a) Controller-Encoder Strategy. The following controller-encoder strategy² achieves the characterization of FTH information CC-Capacity. ( ). ² For any square matrix D with real entries D^(1/2) is its square root. A _(i) ^(g,*) =ē* _(i)(X _(i) ,Y ^(g,*,i-1))=Γ*_(i,i-1) Y _(i-1) ^(g,*)+Δ*_(i) {X _(i) −{circumflex over (X)} _(i|i-1)}, Δ*_(i) =K _(Z) _(i) ^(*,1/2)Σ_(i|i-1) ^(−1/2),Δ*_(i)>0, Y _(i) ^(g,*)=(C _(i,i-1) +D _(i)Γ*_(i,i-1))Y _(i-1) ^(g,*) +D _(i)Δ*_(i) {X _(i) −{circumflex over (X)} _(i|i-1) }+V _(i) ,i=0,1, . . . ,n.

Moreover, from the properties of Kalman filter, the following hold.

(b) Filter Estimates. The innovations process defined by {v*_(i)

Y_(i) ^(g,*)−E{Y_(i) ^(g,*)|Y^(g,*,i-1)}:i=0, . . . , n} satisfies v* _(i) =Y* _(i)−(C _(i,i-1) +D _(i)Γ*_(i,i-1))Y _(i-1) ^(g,*) =D _(i)Δ*_(i) {X _(i) −{circumflex over (X)} _(i|i-1) }+V _(i) ,i=0, . . . ,n. E{v* _(i) |Y ^(g,*,i-1) }=E{v* _(i)}=0,i=0, . . . ,n, E{v* _(i)(v* _(i))^(T) |Y ^(g,*,i-1) }=D _(i) K* _(Z) _(i) D _(i) ^(T) +K _(V) _(i) =E{v* _(i)(v* _(i))^(T)}

and the sequence of RVs, {v*_(i): i=0, . . . , n} is uncorrelated.

The optimal filter estimates and conditional covariances satisfy the following recursions.

$\begin{matrix} \begin{matrix} {\mspace{79mu}{{\hat{X}}_{{i + 1}|i} = {{A_{i}{\hat{X}}_{i|{i - 1}}} + {\Psi_{i|{i - 1}}\left\{ {Y_{i}^{0,*} - {\left( {C_{i,{i - 1}} + {D_{i}\Gamma_{i,{i - 1}}^{*}}} \right)Y_{i - 1}^{\theta,*}}} \right\}}}}} \\ {{= {{A_{i}{\hat{X}}_{i|{i - 1}}} + {\Psi_{i|{i - 1}}v_{i}^{*}}}},{i = 0},\ldots\mspace{14mu},n,{{\hat{X}}_{0|{- 1}} = {given}},} \end{matrix} & \; \\ {{\sum\limits_{{i + 1}|i}{= {{A_{i}{\sum\limits_{i|{i - 1}}A_{i}^{T}}} - {A_{i}{\sum\limits_{i|{i - 1}}{{\left( {D_{i}\Delta_{i}^{*}} \right)^{T}\left\lbrack {{D_{i}K_{Z_{i}}^{*}D_{i}^{T}} + K_{V_{i}}} \right\rbrack}^{- 1} \cdot \left( {{{D_{i}\Delta_{i}^{*}{\sum\limits_{i|{i - 1}}A_{i}^{T}}} + {G_{i}K_{W_{i}}G_{i}^{T}}},\sum\limits_{0|{- 1}}} \right)}}}}}},\mspace{79mu}{\sum\limits_{0|{- 1}}{= {given}}}} & \; \end{matrix}$

where the filter gains {Ψ_(i|i-1): i=0, . . . , n} and the output process {Y_(i) ^(g,*): i=0, . . . , n} are defined by Ψ_(i|i-1)

A_(i)Σ_(i|i-1)(D _(i)Δ*_(i))^(T)[D _(i) K* _(Z) _(i) D _(i) ^(T) +K _(V) _(i) ]⁻¹ Y _(i) ^(g,*)

(C _(i,i-1) +D _(i)Γ*_(i,i-1))Y _(i-1) ^(g,*) +v* _(i) ,i=0, . . . ,n.

Moreover, the a-algebra generated by {Y_(k) ^(g,*): k=0, 1, . . . , i} denoted by

_(0,i) ^(Yg,*)

σ{Y₀ ^(g,*), Y₁ ^(g,*), . . . , Y_(i) ^(g,*)} satisfies

_(0,i) ^(Yg,*)=

_(0,i) ^(v*)

σ{v*₀, v*₁, . . . , v*_(i)}, i=0, . . . , n.

(c) Realization of Optimal Randomize Strategy of L. The controller-encoder strategy {ē*_(i)(⋅,⋅): i=0, 1, . . . , n} realizes the optimal randomized strategy, that is P ^(ē*)(A _(i) ^(g,*) ∈da _(i) |y _(i-1))=π_(i) ^(g,*)(da _(i) |y _(i-1))˜N(Γ*_(i,i-1) Y _(i-1) ^(g,*) ,K* _(Z) _(i) ),i=0, . . . ,n.

(d) FTH Information CC-Capacity Achieving Contoller-Encoder Strategy. The strategy {ē*_(i)(⋅,⋅): i=0, . . . , n} achieves the FTF information CC-Capacity, that is, the following identities hold.

$\begin{matrix} {{\lim\limits_{n\rightarrow}{\frac{1}{n + 1}{I\left( X^{n}\rightarrow Y^{g,8,n} \right)}}} = {\sum\limits_{i = 0}^{n}{I\left( {X_{i};\left. Y_{i}^{g,*} \middle| Y^{g,{*{,{i - 1}}}} \right.} \right)}}} \\ {= {\lim\limits_{n\rightarrow}{\frac{1}{n + 1}{\sum\limits_{i = 0}^{n}{I\left( {A_{i}^{g,*};\left. Y_{i}^{g,*} \middle| Y^{g,{*{,{i - 1}}}} \right.} \right)}}}}} \\ {= {\lim\limits_{n\rightarrow}{\frac{1}{n + 1}{\sum\limits_{i = 0}^{n}\left\{ {{H\left( v_{i}^{*} \right)} - {H\left( V_{i} \right)}} \right\}}}}} \\ {= {{J_{{A\;}^{\infty}\rightarrow Y^{\infty}}\left( {\pi^{g,*},\kappa} \right)} = {C(\kappa)}}} \end{matrix}$

CC-Capacity Achieving Controller-Encoder

P. Consider an example Gaussian-G-RCM-2 of M. The method of O. can be repeated by replacing {(Γ*_(i,i-1), K*_(Z) _(i) )=i=0, 1, . . . , n−1} in (a)-(d) of O. by (Γ^(∞,*), K*_(Z)) obtained in M.

Q. Consider the example of N.

(a) Gaussian RV Message X˜N(0, σ_(X) ²). Then the message is X_(i+1)=X_(i), X₀=X, i=0, . . . , n−1. The Rate Distortion Function (RDF) of the Gaussian RV X subject to Mean-Square Error distortion is given by

${R(\Delta)} = {{\inf\limits_{\hat{X}\overset{\Delta}{=}{{E|X} = {\hat{X}|^{2}{\leq \Delta}}}}{I\left( {X;\hat{X}} \right)}} = {\frac{1}{2}\log\;{\max\left( {1,\frac{\sigma_{X}^{2}}{\Delta}} \right)}}}$

Let {circumflex over (X)}_(i|i)

E{X|Y^(g,*,*)}, i=0, . . . , n be the decoder of message X. Calculating the time-invariant scalar version (i.e., p=q=1) of optimal controller-encoder O., after (n+1) times uses of the control system, the Mean-Square Error (MSE) of the decoder decays geometrically, according to the expression E|X−{circumflex over (X)} _(n|n)|²=σ_(X) ² e ^(−2C) ^(0,n) ^((κ))

where for large enough n, we have C_(0,n)(κ)=(n+1)C(κ), and C(κ) is the capacity of the control system given in N. for cases, |C|>1, |C|<1 (stable and unstable). Letting Δ=Σ_(n|n) we obtain R(Δ)=C _(0,n)(κ).

This implies the controller-encoder-decoder strategy meets the RDF with equality and no other controller-encoder-decoder scheme, no matter how complex it can be, can achieve a smaller MSE. Moreover, lim_(n→∞)E|X−{circumflex over (X)}_(n|n)|²=0. We can also compute the error probability for any finite n.

(b) Equiprobable messages X∈{0, 1, . . . , M_(n)}, n=0, 1, . . . Similar to (a), we can apply the Schalkwijk-Kailath coding scheme, to show for M_(n)

exp{(n+1)R} such that

${R < {\frac{1}{n + 1}{C_{0,n}(\kappa)}}},$ the probability of Maximum Likelihood (ML) decoding error decreases doubly exponentially in (n+1).

(c) The analysis of the multi-dimensional versions is also feasible using the method described.

R. In the method of L.-Q.

(i) the optimal controller-encoder strategy controls the controlled process {Y_(i): i=0, . . . , n} and encodes the information process {X_(i): i=0, . . . , n}, and

(ii) the information process is reconstructed at the output of the control system, using a minimum mean-square error decoder, via the Kalman filter.

(iii)

${C(\kappa)} = {{\lim_{n\rightarrow\infty}{\frac{1}{n + 1}{J_{A^{n}\rightarrow Y^{n}}\left( {\pi^{g,*},\kappa} \right)}}} = {\lim_{n\rightarrow\infty}{\frac{1}{n + 1}{I\left( X^{n}\rightarrow Y^{g,{*{,n}}} \right)}}}}$ is a variant of per unit time infinite horizon stochastic optimal control theory; it reveals that optimal randomized strategies are capable of information transfer, over the control system, which acts as a feedback channel, precisely as in Shannon's operational definition of Capacity of Noisy Channels.

S. The methods described in this patent, for simultaneous control and encoding of information apply to any dynamical system model with inputs and outputs, such as, biological control models, financial portfolio optimization and hedging models, quantum dynamical control systems, etc. 

What is claimed is:
 1. A computer-implemented method for simultaneously controlling a control system and transmitting information signals in the control system, the method comprising: receiving, at a controller of a control system, an information signal to be propagated through the control system, where the information signal includes data to be transmitted through the control system; determining, at the controller of the control system, a control signal to control one or more operations in the control system; encoding, at the controller of the control system, the received information signal to produce an encoded information signal; determining, at the controller of the control system, a two-part signal comprising: 1) the control signal, and 2) the encoded information signal, where the two-part signal is based on a transmission rate below a control-coding capacity of the control system, where the control-coding capacity is a maximum rate of transmitting signals through the control system; transmitting, by the controller, the two-part signal using the control system as a communication channel.
 2. The computer-implemented method of claim 1, wherein encoding the information signal, further includes: encoding parameter values, or improper or failure regimes of actuator devices.
 3. The computer-implemented method of claim 1, wherein encoding the information signal message further includes: removing redundancy in data.
 4. The computer-implemented method of claim 1, further comprising: receiving, at one or more devices of the control system, a signal that conveys information about the information signal; and extracting, at the one or more devices, the information signal.
 5. The computer-implement method of claim 4, further comprising: determining, at one or more devices of the control system, feedback for the controller; and transmitting, by the one or more devices of the control system a process output signal comprising the encoded information signal.
 6. The computer-implemented method of claim 5, further comprising: receiving, at a decoder, the process output signal; and decoding, at the decoder, the encoded information based on an encoder scheme.
 7. The computer-implemented method of claim 6, wherein the two-part message indicates the encoder scheme.
 8. The computer-implemented method of claim 5, further comprising: receiving the process output signal at one or more of: 1) a controller, or 2) a decoder device.
 9. The computer-implemented method of claim 1, wherein the control system is a dynamic system, and further includes a communication channel with memory.
 10. The computer-implemented method of claim 1, wherein encoding the received information signal comprises: determining the control-coding capacity of the control system based on a system model of the control system; and implementing a coding scheme to encode the information signals based on the determined control-coding capacity.
 11. The computer-implemented method of claim 10, wherein the system model of the control system comprises one or more of: power constraints, directed information pay-off, noise distributions, cost functions, feedback models, or an information rate of the control system.
 12. The computer-implemented method of claim 10, wherein determining the control-coding capacity includes solving the optimization of a representation of the control-coding capacity based on the system model, wherein the representation of the control-coding capacity includes one or more conditional control input distributions and one or more cost functions.
 13. The computer-implemented method of claim 10, wherein transmitting the two-part signal further comprises: transmitting the two-part signal based on the control-coding capacity to achieve optimal performance of one or more control objectives.
 14. The computer-implemented method of claim 1, wherein the control system is a dynamic system including one or more inputs, one or more outputs, and one or more controllers.
 15. A system for simultaneously controlling a dynamic system and transmitting information signals in the dynamic system, the system comprising: a controller for generating and transmitting control signals to control one or more operations of a dynamic system; an encoder communicatively coupled to the controller, one or more non-transitory storage media configured to store processor-executable instructions; and one or more processors operatively connected to the controller and the one or more non-transitory storage media and configured to execute the processor-executable instructions to cause the system to: receive, at the controller, an information signal to be propagated through the dynamic system, where the information signal includes data to be transmitted through the dynamic system; determining, at the controller, a control signal to control one or more operations in the dynamic system; encoding, at the encoder, the received information signal to produce an encoded information signal; determining, at the controller, a two-part signal comprising: 1) the control signal, and 2) the encoded information signal, where the two-part signal is based on a transmission rate below a control-coding capacity of the dynamic system, where the control-coding capacity is a maximum rate of transmitting signals through the dynamic system; transmitting, by the controller, the two-part signal using the dynamic system as a communication channel. 